Rabi spectroscopy of a qubit-fluctuator system 
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Superconducting qubits often show signatures of coherent coupling to microscopic two-level fluc- 
tuators (TLFs), which manifest themselves as avoided level crossings in spectroscopic data. In this 
work we study a phase qubit, in which we induce Rabi oscillations by resonant microwave driving. 
When the qubit is tuned close to the resonance with an individual TLF and the Rabi driving is 
strong enough (Rabi frequency of order of the qubit-TLF coupling), interesting 4- level dynamics 
are observed. The experimental data shows a clear asymmetry between biasing the qubit above or 
below the fluctuator's level-splitting. Theoretical analysis indicates that this asymmetry is due to 
an effective coupling of the TLF to the external microwave field induced by the higher qubit levels. 
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Spectroscopic analysis of superconducting qubits often 
shows clear signatures of avoided level crossings, indi- 
cating the presence of microscopic two-level fluctuators 
(TLFs) that can be in resonance with the qubit. Evi- 
dence for the existence of TLFs have been found in nearly 
all known types of superconducting qubits, including 
phase-Cl^, fiux-l^, charge-!^, and transmon qubit^. Since 
TLFs are considered to be a source of decoherencd^^^, 
experiments are usually conducted by biasing the qubit 
in a frequency range where none of these strongly coupled 
natural two- level systems are present. Alternatively, one 
can take advantage of the longer coherence times of TLFs 
as compared to the qubits for using them as a quantum 
memory^. Here we focus on the dynamics of the qubit- 
fluctuator system on or near resonance. 

There are at least two possible mechanisms explaining 
the interaction of the TLFs with the qubit: (i) the TLF 
is an electric dipole which coup les t o the electric field 
in the qubit 's Josephson junctioif^. Nanoscale dipoles 
could emerge from metastable lattice configurations in 
the amorphous dielectric of the junction's tunnel bar- 
rieilH^ (ii) the state of the TLF affects the critical cur- 
rent of the qubit's Josephson junction^ In this case the 
TLF could be related, e.g., to the formation of Andreev 
bound states at t he int erface between the superconductor 
and the insulatoi^^. 

In this Letter, we explore the complexity of the dy- 
namical behavior of a driven phase qubit operated in the 
vicinity of a resonance with a two-level fluctuator. Due 
to the strong coupling between the qubit and the TLF 
and equally strong Rabi driving, we observe the dynam- 
ics of the resulting 4-level hybrid system consisting of 
the microscopic defect state and the macroscopic phase 
qubit. Strong microwave driving of the coupled system 
leads to coherent oscillations, revealing a characteristic 



beating pattern which we analyze quantitatively. Our 
experimental data displays a distinct asymmetry of the 
system response with respect to the detuning between 
the qubit and the TLF. We argue that this asymmetry 
is due to Raman-like virtual processes involving higher 
quantum levels of the qubit, giving rise to an effective 
driving of the TLF. 

The sample investigated in this study is a phase qubilPl, 
consisting of a capacitively shunted Josephson junction 
embedded in a superconducting loop. Its potential en- 
ergy has the form of a double well for suitable combina- 
tions of the junction's critical current (here, Ic = 1.1/iA) 
and loop inductance (here, L = 720 pH). For the qubit 
states, one uses the two Josephson phase eigenstates of 
lowest energy which are localized in the shallower of the 
two potential wells, whose depth is controlled by the ex- 
ternal magnetic flux through the qubit loop. The qubit 
state is controlled by an externally applied microwave 
pulse, which in our sample is coupled capacitively to the 
Josephson junction. A schematic of the complete qubit 
circuit is depicted in Fig.[TJa). Details of the experimen- 
tal setup can be found in Ref. [Ml During all measure- 
ments presented in this paper, the sample was cooled to 
a temperature of 35 mK in a dilution refrigerator. 

Spectroscopic data taken in the whole accessible fre- 
quency range between 5.8 GHz and 8.1 GHz showed only 
4 avoided level crossings due to TLFs having a coupling 
strength larger than 10 MHz, which is about the spec- 
troscopic resolution given by the linewidth of the qubit 
transition. In this work, we studied the qubit interact- 
ing with a fluctuator whose energy splitting was Cf /h = 
7.805 GHz. From its spectroscopic signature shown in 
Fig.jljb), we extract a coupling strength v±/h ^ 25 MHz. 
The coherence times of this TLF were measured by di- 
rectly driving it at its resonance frequency while the qubit 
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was kept detunecF^. A tt pulse was applied to measure 
the energy relaxation time Tij ^ 850 ns, while two de- 
layed 7r/2 pulses were used to measure the dephasing time 
T2 f ~ 110 ns in a Ramsey experiment. To read out the 
resulting TLF state, the qubit was tuned into resonance 
with the TLF to realize an iSWAP gate, followed by a 
measurement of the qubit's excited state. 
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FIG. 1: (color online) (a) Schematic of the phase qubit cir- 
cuit, (b) Probability to measure the excited qubit state 
(color-coded) vs. bias flux and microwave frequency, reveal- 
ing the coupling to a two-level defect state having a resonance 
frequency of 7.805 GHz (indicated by a dashed line). 

Experimentally, we observe the probability P{\e)) of 
the qubit being in its excited state after driving it res- 
onantly with a short microwave pulse. Varying the du- 
ration r of the microwave pulse allows us to observe the 
evolution of P{\e)) in the time domain. If the energy 
splitting of the qubit is tuned far away from that of the 
fiuctuator, the qubit remains decoupled from the TLF 
and P{\e)) displays the usual Rabi oscillations in the form 
of an exponentially decaying sinusoid having only a sin- 
gle frequency component. For our qubit sample, which 
has coherence times of Ti^q ^ 120 ns and T2^q ~ 90 ns, 
these oscillations have the characteristic decay time of 
about 115 ns. If, in contrast, the qubit is tuned close 
to the resonance frequency of a TLF, the probability to 
measure the excited qubit state shows a complicated time 
dependence, which is very sensitive to the chosen qubit 
bias. 

Figure |2|a) shows a set of time traces of P{\e)) taken 
for different microwave drive frequencies. Each trace was 
recorded after adjusting the qubit bias to result in an 
energy splitting resonant to the chosen microwave fre- 
quency. The Fourier transform of this data, shown in 
Fig.[2|b), allows us to distinguish several frequency com- 
ponents. Frequency and visibility of each component de- 
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FIG. 2: (color online) (a) Experimentally observed time evo- 
lution of the probability to measure the qubit in the excited 
state, P(|e))(t), vs. driving frequency; (b) Fourier-transform 
of the experimentally observed P(|e))(t). The resonance fre- 
quency of the TLF is indicated by vertical lines, (c) Time 
evolution of P{\e)) and (d) its Fourier-transform obtained by 
the numerical solution of Eq. ([sj as described in the text, tak- 
ing into account also the next higher level in the qubit. (As 
the anharmonicity A/h ^ 100 MHz in our circuit is relatively 
small, this required going beyond the second order perturba- 
tion theory and analyze the 6- level dynamics explicitly). The 
qubit's Rabi frequency Qq/h is set to 48 MHz. 



pend on the detuning between qubit and TLF. We note a 
striking asymmetry between the Fourier components ap- 
pearing for positive and negative detuning of the qubit 
relative to the TLF's resonance frequency, which is indi- 
cated in Figs. |2|a,b) by the vertical lines at 7.805 GHz. 
We argue below that this asymmetry is due to virtual 
Raman-transitions involving higher levels in the qubit. 
To describe the system theoretically, we write down 
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the Hamiltonian, consisting of two parts: H = Hs + 
Hi, with Hs being the system Hamiltonian, representing 
qubit, TLF and their couphng and Hj describing the 
interaction between system and microwave driving. The 
Hamiltonian of the qubit circuit reads 



Ecin-nc) - Ej cos cj) ^ El {(j) - 



t)\ (1) 



where Ec/j/l are charging/ Josephson/inductive energies 
of the circuit, (j) is the phase difference across the Joseph- 
son junction, and n is the dimensionless charge conjugate 
to (j), i.e., [0, n] = i. The circuit can be manipulated by 
applying an ac driving to gate charge no or the external 
flux (pext- The TLF is described as a two level system 
Hg = l/2efTz which couples either to the electric field 
across the junction ex {n — no) or, alternatively, to the 
Josephson energy ex cos^. The coupling can be either 
transverse, oc t±, or longitudinal, cx r^. 

For maximum generality, we first define a minimal 
model needed to describe the splitting of Fig. [l] To this 
end, we restrict ourselves to the lowest two states of the 
phase qubit circuit (the qubit subspace) and disregard 
the longitudinal coupling cx . Within the rotating wave 
approximation (RWA) the Hamiltonian reads 



Tjmin 



^q^z + + + cr+r_) , (2) 



with the Pauli-matrices for the qubit a and for the fluc- 
tuator r. The minimal interaction Hamiltonian couples 
only the qubit to the driving field via the coupling con- 
stant Qqi Hf'"" = Qq COS (udt) (7^. The RWA is justified 
since ^ v± ^ ^ ^f- Rabi oscillations in this mini- 
mal system have been considered earliei . 

Going to the rotating frame for both qubit and TLF 
and taking the frequency of the driving to be resonant 
with the qubit splitting, uj^ = we arrive at the effective 
4-level Hamiltonian 

J?™" = -^SajT, + ^vi_ (a_r+ + (t+t_) + ^n^a, , (3) 

where Suj = {cq — Cf). The level structure and the spec- 
trum of possible transitions in the Hamiltonian (|3| is 
illustrated in Fig. [3^. The transition frequencies in the 
rotating frame correspond to the frequencies of the Rabi 
oscillations observed experimentally. 

To describe the time evolution of our system we con- 
sider the state = ^kCke~'^^^^ \k) , where Ek are 
the eigenvalues and \k) the eigenstates of the Hamilto- 
nian ([3]). The coefficients Ck are determined by the initial 
conditions. The eigenvectors \k) can be expressed as lin- 
ear combinations of the measurement basis states \g \), 
\g^), let), i-e., the mutual eigenstates of and 
r^, which we denote by {|/)} with / = 0, 1,2,3. For the 
expectation value of the operator we get 



(^k.lO , (4) 



where a^^i = Ck {I \k) and we used the fact that cr^ is di- 
agonal in basis {|/)}- From Eq. Q we can extract the 
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FIG. 3: (color online) (a) Analytically obtained transition 
spectrum of the Hamiltonian (|3| in the minimal model for 
= 40 MHz and v^jh = 25 MHz. Dashed-dotted lines 
show the transition frequencies while the gray-scale inten- 
sity of the thicker lines indicates the weight of the respec- 
tive Fourier-components in the probability P(|e)). The sys- 
tem shows a symmetric response as a function of the detun- 
ing buj. Two of the four lines are double degenerate, (b) 
The same as (a) but including the second order Raman pro- 
cess with = v\_Vlql IS.. The two degenerate transitions in 
(a) split and the symmetry of the response is broken. In- 
set: Schematic representation of the structure of the Hamil- 
tonian (|3|. We denote the ground and excited states of the 
qubit as"^') and |e) and those of the TLF as ||) and |t)- Ar- 
rows indicate the couplings between qubit and fluctuator 
and to the microwave field Qg and Qj. 



Fourier coraponents of the experinientally raeasured ex- 
cited state population P(|e)) = (1 + {(Jz))/2. There are 
six coDiponents with, in general different, transition fre- 
quencies Era — Ek. These are shown in Fig. [3^ for the 
minimal model. Only four lines are seen due to two dou- 
ble degeneracies. The intensity of the thick lines overlay- 
ing the dashed-dotted transition lines corresponds to the 
amplitude of these Fourier components. The situation 
depicted in Fig. [3] and realized in our experiment corre- 
sponds to the qubit and the fluctuator initially in their 
ground states. It is important to note that the pattern 
of Fig. |3]is characteristic for the regime Vtq v^. 

As seen in Fig. [3^ the observed asymmetry in the re- 
sponse can not be explained by the minimal model. We 
identify three possible mechanisms which could break the 
symmetry: (i) Longitudinal coupling between qubit and 



TLF H 



long 



v\\(JzTz. We note that the longitudinal cou- 



pling is excluded for the electric dipole coupling mech- 
anism in phase and flux qubits, since this term would 
necessitate an average electric field (voltage) across the 
junction. The longitudinal coupling might be present if 
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the TLF couples via a change in the critical currentPE^. 
In this case the state of the TLF directly affects the 
shape of the Josephson potential, therefore modulating 
the level-splitting of the qubit. For realistic parame- 
ters, this might lead to a strong longitudinal coupling 
v\\ . Such a coupling was, however, ruled out spectroscop- 
ically in Ref. 4 as well as by our preliminary spectroscopic 
dat#^. (ii) Direct coupling of the TLF to the external 
field Hj = cos (udt) r^- Due to the presumably small 
size of the TLF this coupling should be negligible, (iii) 
Effective coupling of the TLF to the external driving field 
due to a second order Raman-like process in which the 
next higher level of the qubit |e2) is virtually excited fol- 
lowed by a mutual flip of the TLF and the qubit (back 
to state |e)). The energy difference between the states 
1 62) and |e) is given by — A, where A characterizes 
the anharmonicity of the qubit. This gives an effective 
coupling Hj = Q'j cos {uJdt) |e) (e|, i.e., the coupling is 
present only when the qubit is excited. For (^cj < A, we 
find ^ v±_rtq//S.. In Fig |3]3 we show the spectrum 
of transitions with only the term added to the min- 
imal model, Eq. ([3|, not including longitudinal coupling 
or direct coupling of the TLF to the driving field. 

To fully describe the experiment, we include decoher- 
ence in our calculations by solving the time evolution of 
the system's density matrix p using a standard Lindblad- 
approach^^. The dynamic equations are given by 



+ U,pL] - \ {l,l\,P^ , (5) 



where the sum is over all possible channels of decoherence 



with the respective rates Fj. The are the operators 
corresponding to each decoherence channel, e.g., pure de- 
phasing of the qubit is described by the operator . The 
theoretical spectral response of the system obtained by 
numerically solving the dynamical equations is shown in 
Fig. [2|c, d). Relaxation and pure dephasing rates for 
qubit and TLF have been taken to be equal to the values 
mentioned earlier. The plot of Fig. |2|c,d) takes into ac- 
count the third level in the qubit. As the anharmonicity 
/S./h 100 MHz is known from other measurements^, 
we have no additional fit parameters and quantitatively 
reproduce the experimental data. Note that we are able 
to explain the experimental data by assuming v\\ = 0, 
which provides further evidence in favor of the dipole 
coupling mechanism. 

In conclusion, we studied the dynamics of a driven sys- 
tem consisting of a phase qubit strongly coupled to a 
TLF. The Fourier- analysis of the Rabi oscillation data 
reveals the characteristic pattern of transition frequen- 
cies in the coupled system. This asymmetric pattern is 
reproduced quantitatively by the presented theoretical 
model including virtual transitions to the qubit 's higher 
levels. The apparent absence of the longitudinal coupling 
between the qubit and the TLF gives a hint about the 
microscopic nature of the TLFs. 
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